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Abstract 

The possibility of a gravitational phase transition, especially with re- 
spect to neutron stars is investigated. First, a semiclassical treatment 
is given, predicting a gravitational London penetration depth of 12km 
for neutron stars. Second, the problem is considered from a Ginzburg- 
Landau point of view. A gravitational Meissner effect, a gravitational 
Aharanov-Bohm type effect and a gravitational ferromagnetic type phase 
are predicted. Finally, a field theoretic consideration predicts a mass term 
for the graviton below a certain critical temperature. 

PACS number(s): 04.60, 04.70, 11.15.E, 97.60.J 
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1 Introduction 



In the post-Newtonian approximation, gravity looks like an effective U(l) theory, 
closely resembling the theory of electrodynamics. In this work we consider 
effects arising from the breaking of this U(l) symmetry. Basically we consider 
the possibility of a macroscopic quantum gravity effect in analogy to the effect 
of superconductivity in the electromagnetic theory. 

Using a semiclassical approach in close analogy to London's original ap- 
proach, we find that for neutron stars, macroscopic quantum gravity effects 
seem to become significant. We predict a gravitational Meissner effect with a 
gravitational London penetration depth for neutron stars of about 12km, which 
is somewhat smaller than the neutron star's radius. 

A Ginzburg-Landau treatment leads to the same result. We use the close 
analogy between electromagnctism and the post-Newtonian approximation and 
predict again the Meissner effect, a gravitational Aharanov-Bohm effect and a 
gravitational ferromagnetic type phase, where all the spin angular momenta of 
the neutron star become aligned. 

Finally we look at the graviton. We use the gravitational wave expansion 
from general relativity, treating the graviton as a spin two field. Adding the 
Ginzburg-Landau terms through the energy momentum tensor, we observe that 
the theory undergoes a phase transition, with the graviton becoming massive 
below a certain critical temperature. Again, this predicts a Meissner effect, 
justifying the semi-classical treatment. 

2 Motivation 

Neutron stars are extraordinarily fascinating objects: Not only are they our 
preferred laboratory to test general relativity, but they also seem to exhibit 
macroscopic quantum effects, i.e., it is believed that they are superfluid and 
maybe superconducting p4j . Although these two issues seem to be separate, we 
will consider here the possibility of a macroscopic gravitational quantum effect, 
in close analogy to the case of superconductivity. 

Neutron stars are made mostly out of neutrons. They have usual radii of 
about 15 km and masses of about 1.4 times the mass of the sun. Theoretical 
calculations indicate that for low enough temperatures (< 10 9 K) the neutrons 
are expected to pair up | po| , similar to nucleons forming pairs inside nuclei in 
the Interacting Boson Model (IBM) Jl^, [2l{| . The calculations indicate that 
in the crust region the singlet 1 5'o state should be the preferred one, however, 
for the majority the triplet z Pi state should be the preferred state p], pj, |32fl. 
These neutron pairs are bosons by statistics and bose condensation is expected. 
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Interesting enough, there is experimental evidence for the pairing: one, it comes 
from the theory of cooling of neutron stars, where the reduction of the specific 
heat due to extensive pairing is essential p3| , and two, evidence for superfluidity 
can be found in the glitches in the timing history of pulsars |^4|, [3(J . 

Therefore, we may assume that there exists a bose condensate inside a neu- 
tron star. This bose condensation then is responsible for the breaking of the 
translational symmetry leading to superfluidity [fyl and possibly the electro- 
magnetic U(l) symmetry leading to superconductivity [jlt^ [3C||. Since gravity in 
the post-Newtonian approximation is an effective U(f ) theory, we may wonder 
what would happen if in this approximation, this gravitational U(l) symmetry 
were to be broken. 



3 Semiclassical Treatment 

For weak gravitational fields and low velocities, Einstein's field equations can 
be written in the post-Newtonian approximation |20[ |37). We are mostly in- 
terested in magnetic-type gravity, and therefore we shall use the truncated and 
rewritten version of the parametrized-post-Newtonian (PPN) formalism given 
by Braginsky et. al. [Q. In this formalism, Einstein's field equations can be 
rewritten in a form very much resembling that of Maxwell's equations 

1 r)TT 

V • g = -47TPQ + O (c- 2 ) , Vxg = --— (1) 

V.H = 0, VxH = 4(-4 7 r P0 I + i|) ) (2) 

where po is the density of rest mass in the local rest frame of the matter and v is 
the ordinary velocity of the rest mass relative to the PPN coordinate frame. The 
equations of motion of an uncharged particle are identical to the electromagnetic 
Lorentz force law ||, 

^=g + ivxH + 0( c - 2 ). (3) 

In this weak field, slow motion expansion of general relativity, g contains 
mostly first order corrections to flat space-time, and H contains second order 
corrections. For the solar system, g is just the normal Newtonian gravitational 
acceleration, whereas H is related to angular momentum interactions and effects 
due to H are about 10 12 times smaller than those due to g 

With these gravitational Maxwell equations we can derive the equivalent 
of London's equations JTJ, |l5|, for gravity. The derivation is completely 
analogous to the electromagnetic case [^3|. This leads to the second London 
equation for gravity 
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Vxj = -VG^H, (4) 

which in turn predicts the Meissner-Ochsenfeld effect |ll| with a gravitational 
London penetration depth of 

A * -(ass) V2 = 5 ' 18x 1012 (5) 

where the density po is measured in kg/m 3 and is given in meters. 

Calculating the London penetration depth for neutron stars with approxi- 
mate densities of about p^s = 2x 10 17 kg/m 3 , we obtain a London penetration 
depth of 12km, which is slightly smaller than the radius of a neutron star with 
that density. The London length is inversely proportional to the density which 
means increasing density leads to decreasing London penetration depth. 

As a result, the onset of a gravitational Meissner effect is predicted. In 
essence it implies that the gravitational magnetic field caused by the huge an- 
gular momentum of the star, will be expelled from the center of the neutron 
star. This could be accomplished through induced matter-supercurrents in the 
outer layers of the neutron stars, creating counter-magnetic fields to expel the 
gravitational magnetic field from its interior, in analogy to the electromagnetic 
case. The gravitational Meissner effect may lead to significant modifications 
with respect to the collapse of a neutron star. 

For our derivation of the London equations to be valid we must require that 
the London penetration depth A^, is much larger than the coherence length £ 
Q. A simple upper estimate for the coherence length yields £o < 10~ 13 m p3| , 
which is larger than the separation between two nucleons, but much smaller 
than the London penetration depth. 

4 Gravitational Ginzburg-Landau Theory 

The previous approach was basically classical in nature. For a Ginzburg-Landau 
treatment fll|l of the problem we need a quantum theory. We first show how 
to modify Schrodinger's equation for the post-Newtonian approximation. Then 
we describe the neutron pair bose condensate. With this in place a standard 
Ginzburg-Landau treatment leads to the prediction of Meissner effect and a 
Aharanov-Bohm type effect. In addition for the triplet state we predict a macro- 
scopic ordering of the spins similar to ferromagnetism. 

4.1 Gravitational Schrodinger equation 

The first experiment showing the direct influence of gravity on the quantum 
mechanical phase of a wave function was done a little more than twenty years 
ago by Colella, Overhauser and Werner, now known as the COW-experimcnt 
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[g| ^5|. They used thermal neutrons in an interferometer type experiment and 
showed that rotating the experimental apparatus in the earth's gravitational 
field gives an interference pattern exactly as predicted by the 'gravitational' 
Schrodinger equation (see Sakurai p8fl), 

-|-V 2 +m^V = ^. (6) 



2m J dt 

Since then additional experiments have been done, as for example the influ- 
ence of the earth's rotation on the phase (Sagnac effect) has been verified ex- 
perimentally |26| . Recently atomic interferometry has become available, which 
allows for much higher precision ||. In the near future this should allow for 
the measurement of such effects as coupling of a particle's spin to the earth's 
rotation, as well as second order correction to the Newtonian potential caused 
by curvature, and the somewhat more distant future, Lense-Thirring effect and 
gravitational waves. 

Although these experiments have not been done yet, we are encouraged to 
follow the analogy of the electromagnetic treatment and propose the Schrodinger 
equation for post-Newtonian gravity, 



where we identify the minimal coupling 



V->V-^A. (8) 

nc 

Here, q = \J~Gm, and details for the derivation and units can be found in the 
appendix. Please note, that equation ([?]) is invariant under the following gauge 
transformation: 

A^A + VA(ir) and ip e iqA W/ Bc V>. (9) 

One major difference to electrodynamics is that (cj>, A) is not a four-vector, 
but rather transforms as a particular component of a rank two tensor, which 
becomes clear considering its definition in terms of the metric tensor, 

.goo = l-20, gij = (1 - 2(f))Sij, and g 0i = -A ; . (fO) 

However, in the following construction we never require it to be a four vector. 



4.2 Long Range Order 

As indicated in the introduction, there is experimental evidence that the neu- 
trons inside a neutron star are superfluid. This means that there exists a bose- 
condensate already, which can be described by a Ginzburg-Landau type theory. 
Let us review the basic ideas behind this long-range ordering. 
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Neutrons are fermions with spin of 1/2. Due to the Pauli exclusion principle, 
no one-particle level can have a macroscopic occupation number in a fermion 
system, so long range order in the one-particle density matrix is forbidden. But 
there can be long range order in the two-particle density matrix ]18| : 



risi,r 2 s 2 



P (2) 



(^>i)<(r 2 )*4(r^ s <(rl)>. (11) 



In the ordered phase, when the pair of points ri, r 2 is far removed from the pair 
r[ , r 2 , this becomes 

(r 1 s 1 ,r 2 s 2 p {2) r[s[,r 2 s' 2 j ~ $* iS2 (ri,r 2 ) $ s > iS > 2 (r[,r' 2 ), (12) 

and such ordering is called pairing. 

In general, in a paired system the order parameter $ SlS2 (ri, r 2 ) has a richer 
structure: it is a 2 x 2 matrix in spin space, each element of which is a complex 
function of the difference variable r — T\ — r 2 . Because the fermion field opera- 
tors anti-commute it follows that the order parameter behaves like a two-particle 
wave function under interchange of particles & S is 2 ( r li r 2) — —®s 2 s 1 ( r 2, r i)- It 
also follows that the order parameter transforms like a two-particle wave func- 
tion under rotations in spin space, rotations in position space, or Galilean trans- 
formations ps| . 

The order parameter can be expanded in a basis of four spin states consisting 
of the usual singlet state |0), where the spins of the two neutrons effectively 
cancel, and three triplet states {x^), where the spins add up to one. The singlet 
state is invariant under spin rotations and the triplet states transform like the 
x-, y-, and z-components of a vector. We now can write the order parameter or 
two-particle wave function as 

$ SlS2 (ri,r 2 ) = (si,s 2 \$(r u r 2 )) , (13) 

where 

*(ri,ra) = <Mri,r 2 ) |0) + tp„(r 1 ,r 2 ) M (14) 
is a combination of singlet cj)(ri,r 2 ) and triplet - 0p( r ii r 2) order parameters. 



4.3 Singlet 1 S State 

In the singlet state, the quantum numbers S and L are zero, meaning that the 
order parameter $> SlS2 {ri,r 2 ) = 4>(ri,r 2 ) is a spin scalar and depends only on 
the magnitude of r. This is equivalent to the type of ordering in helium-4, and 
the only degree of freedom of the order parameter is associated with the phase 
degeneracy. 

Since only the two-particle wave function (jxp* has a non-vanishing expecta- 
tion value, physically measurable properties of a spatially uniform system can 
depend only on the product If, in the spirit of Landau, we expand the free 
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energy of ordering in powers of the order parameter, the first two non- vanishing 
terms will give a free energy density of the form 

^o = M 2 |0| 2 + A|0| 4 . (15) 

Slow spatial variations in the order parameter can be treated in the usual 
gradient expansion. To leading order there is only one term invariant under 
rotations and gauge transformations and Tq becomes 

T= l -K\v<t>\ 2 + ^\4>\ 2 + \\4>\ 4 . (16) 

The coupling to gravity can be accomplished through minimal coupling. 
Recall equation (||) and the fact that the mass of the neutron pair is 2q, we find 
for the free energy 

T = f ( v " l Yc A ) * ( v + l ¥c A ) r + + A W * )2 ' (17) 

From here we can obtain the currents associated with the gauge symmetry (j^j 
through variation, remembering that the mass of the neutron pair is 2q. Thus 
we find for the conserved current density 

j = ^(#*V^-#V0*-^A|0| 2 V (18) 

Now assume that <f> = \<p\ e %v has its spatial variation in the phase <p and not 
in the magnitude \<f>\, then the current simplifies to 

Since the curl of any gradient vanishes, and since |</>| 2 is essentially constant, we 
immediately deduce the second London equation 

V x j = — ^ — \4>\ V x A. (20) 
h c 

Comparing this with our result from the previous chapter, equation (^), we 
can use this equation to relate macroscopic to microscopic quantities, that is, 
we can identify the density 

= (21) 



Notice that by redefining <j> — * vK<j) in equation (|16|), we can absorb the K into 
the \4>\ ■ Simultaneously we will also redefine fi — * \fK^i and A — > i^T 2 A. This 
then allows us to express |^>| 2 in terms of known quantities, i.e., 

-2 



H 2 = ^Po- (22) 
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4.3.1 Meissner Effect 



One peculiar feature of the Ginzburg-Landau free energy is that its minimum 
depends on temperature. For temperatures T above a certain critical temper- 
ature T c , the minimum of the potential will simply be at |0| 2 = 0. When the 
temperature T is below the critical temperature T c , i.e., T < T c , then the min- 
imum is at 

|0| 2 = > 0. (23) 

Furthermore, we assume that 4> varies only very slightly over the sample, i.e., 
|0| is essentially constant, 



and taking the curl of Ampere's equation 



Vxj = -dr W Vx A, (24) 
n c 



wc find 



r, x t 167T. 

V x H = j 

c 



\ II rr ( g ) |0| 2 H = fc 2 H, (25) 



where we defined 

For a neutron star of density 2 x 10 17 kg/m 3 we can calculate k using the fact 
that k = 1/A L = 7.4 x 10" 9 l/m 2 . 

4.3.2 Flux Quantization: Gravitational Aharanov-Bohm Effect 

As a result of the Meissner effect, appreciable currents can only flow near the 
surface of a gravitational superconductor. As a consequence, if we integrate the 
current over a closed ring inside the superconductor, we find that 

O = ^j.dl = -||0| 2 ^V^+^A^) -dl. (27) 

Then Stokes theorem gives us 

Adl= J VxAdS = j H • dS = (28) 

where $>h is the flux enclosed in the ring. Since the order parameter $ is 
single- valued, its phase must be an integer multiple of 27r, and therefore 

j)V(p- dl = 27m = = Y c j A dl ( 29 ) 
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or 

he 

<5> H = —n = $ n (30) 
iq 



with the quantity <I>o = hc/2q = 59.8 [m 3 /s 2 vG] being the fluxoid or flux quan- 
tum. 

4.4 Triplet 3 F 2 State 

Only neutrons in the outer layers of a neutron star will be in the singlet state. 
Numerical simulations suggest that most of the neutrons will be found in the 
triplet, the 3 i-2 state ||, ||, |32|. Therefore, we will briefly consider the 3 P2 state 
with quantum numbers S = 1, L = 1 and J = 2. In the triplet state, <j) vanishes 
in equation ( |l4|) and the vector ip^ is of the form 

Mi") =A fiV r vX (\r\), (31) 

where is a complex 3x3 matrix. Under rotations of the orbital coordinates 
A^y transforms as a vector with respect to index is, under rotations of the spin 
coordinates, A^ v transforms as a vector with respect to the index fi. In general, 
A^ represents all Z = 1 order parameters. To stay in the 3 P2 subspace we must 
restrict the matrix A to be traceless and symmetric. Then the most general 3 P2 
Ginzburg-Landau functional to order four, is |22|, |3(], |3^| 

T Q = -a TrAA f + ft | Tr A 2 1 2 + fa (TrAA*) 2 + fa (TrA 2 A* 2 ) . (32) 
o 

This equation represents the most general mean field theory of 3 i-2 pairing near 
T c consistent with the overall rotational and gauge invariance. 

The gradient term for the triplet state looks somewhat more complicated, 



T grad = ^ (v,a;,) (v^ aii )+ : y (v,a; (v,a aii )+ : y (v^;.) ( V ^ M -) 

" (33) 

but it gives a current similar to (|l|), predicting a Meissner effect. Notice, that 
for example for helium-3 the three coefficients are nearly the same K\ = K2 = 
Kz @. 

What changes significantly in the triplet state are the minima of the free 
energy. Again, above a certain critical temperature, the minimum will be for 
J-q = 0. However, below this temperature there are several different configura- 
tions. 

Stability requires that the fourth-order term in the Ginzburg-Landau func- 
tional be positive, which implies that fa > 0, and restricts the fa — fa plane to 
the region fa > ~3fa and fa + 2 fa > —2 fa. Minimizing the free energy leads 
to the following three subregions |ll], |2^] : 



9 



Region I: 03 > \0i\— 0i. In this region the order parameter is unique, except 
for a constant phase factor and a rotation. At the minimum, the free energy is 
given through 

_ a*_ 
min " 4/V 

Region II: > 03 > — 60%. Except for a constant phase factor and a rotation 
the solution is again unique and at the minimum 

J mm 



4 (02 + |/%) 



Region III: 03 < —40% — 2 \0i \ . Here the minimum is given by any A that is 
real, traceless and symmetric, and at the minimum 

J~ min 



4 (01 +02 + §&) 



4.4.1 Macroscopic Quantum Coherence 

While the expectation value of the orbital angular momentum along any axis 
in region II and III vanishes, for the solution in region I it will be maximally 
aligned ]l7| |. This will lead to an effect analogous to ferromagnetism or the 
alignment of spins in helium-3 |55|. Assuming all the orbital angular momenta 
of the 3 P2 state are aligned, the neutron star will have a total, macroscopic spin 
angular momentum of 

S = N-h~^--h = 1.8-10 23 Js. (34) 
m N 

Comparing this to the approximate S = 2 ■ 10 40 Js of the Crab pulsar due to 
rotation, we find that the pulsar's period would have to slow to a period of a 
few hundred million years, for this contribution to become significant. 

Furthermore, one would expect gravitational infinite conductivity to occur. 
For example Weinberg shows that no energy gap is necessary, since infinite 
conductivity depends only on the spontaneous breakdown of gauge invariance 




5 The Graviton becomes Massive 

The graviton is a spin two field. This may be the main objection to the previous 
section. Here we derive the Meissner effect for a spin two graviton. We start 
with the Lagrangian for a scalar field and add a Ginzburg-Landau potential. 
After we have the Lagrangian identified, we proceed to calculate the energy- 
momentum tensor. Expanding the gravitational field in the gravitational wave 
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expansion, we find that the graviton becomes massive below a certain critical 
temperature, i.e., a phase transition has occurred and we expect the Meissner 
effect as a consequence. We then repeat the calculations for a vector field instead 
of a scalar, and arrive at very similar results. We end with a prediction of the 
graviton's mass. 



5.1 Scalar Field 

Let us consider the Lagrangian of a complex scalar field 0, with terms at most 
quadratic in the fields, (and we will leave tadpole like terms out of considera- 
tion) . A good choice is found in Veltman J34| or Feynman j[o) , and also Callan 
et. al. §: 

C = V9 (^g^d^dvp -\(V 2 ~ 2ai?) 00* - A (00*) 2 + bR^d^dvft 

(35) 

We will not consider the b term and also higher terms, since it will imply more 
than two derivatives. In addition, we can redefine /i 2 — > /i 2 — 2aR to simplify 
the above Lagrangian. We find, 

C = V9 f-i^a W - ~ A (00* ) 2 ) ■ (36) 

Through variation with respect to the metric we will obtain the energy- 
momentum tensor for the above Lagrangian, i.e., 

r — (37) 



and we find that the energy momentum tensor is given through 
1 



S k 09a0* - g K \ Q 5 a/3 5 Q 09 (3 0* - i^ 2 00* - A (00*)' 



(38) 



We are especially interested in the weak field approximation used for grav- 
itational waves, i.e., — rf- v — kN 1v . Using this expansion, we find for the 
energy momentum tensor 

T KX = ^ K 09 A 0* (39) 

+\ H«a v af3 + « (vnx h aP - h KX n *) + n 2 h KX h a P) a Q 0^0* 

+\ V.x Qa* 2 #* + A (00*) 2 ) + Q^ 2 00* + A (00*) 2 

Notice that for both possible minimal values of |0| 2 , that is, |0| 2 = above, 
and |0| = — /i 2 /2A below the transition temperature, the vacuum expectation 
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value of the last two terms vanishes. Next for simplicity neglect terms in k 2 and 
furthermore assume that all the space-time dependance is in the phase, i.e., 

d K (f> = d K (|0| e*) = i<f>d K <p. (40) 

We then obtain the following form for the energy momentum tensor 

-d K (pd\Lp + ~~t/ k \ d a ipd a ip - ^-t] k \ h al3 d a Lpdf]Lp - ^-h KX d a ipd a ip 



T K x - 



101* ■ 



(41) 

This gives the left hand side of Einstein's equations, 

R a p + ^r)apR=-^GT aP . (42) 

For the right hand side we use the Ricci tensor in the gravitational wave expan- 
sion m 

R a p = {d x d x h aP + d a dfih\) , (43) 

where we made use of the gauge condition d\h\ — 0. With this we find the 
following wave equation for the graviton 

d a d a h KX + (d K d x + nnxd^) h a a = (44) 



-8ttGk 



K 1 

-d K ipdxV - -t}k\ h a0 d a ipdpip + - (r] KX - nh KX ) d a (pd a (p 



\0\ 
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Now the vacuum expectation value of |(/>| 2 above and below phase transition 
is given by 

|0| 2 = for T>T C (45) 

I0I 2 = for T<T C , (46) 

and we see that above the critical temperature, the graviton is massless, but 
below Tc the graviton has acquired a mass. 

Clearly this is a case of spontaneous symmetry breaking. Analogous to 
the Higgs mechanism we expect the graviton to acquire three new degrees of 
freedom, since it is a spin two field. Those degrees can only come from the scalar 
field, which means we need a scalar triplet or three scalars to conspire together 
to give one graviton its mass. This seems somewhat awkward and artificial, but 
still could be done. However, if the symmetry breaking field is a vector field 
instead of a scalar, this scenario becomes much more natural: One vector field 
has three degrees of freedom, and hence one graviton would eat up one vector 
field to acquire its mass, instead of three scalar fields. Thus we will consider the 
vector field next. 
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5.2 Vector Field 

We have two reasons to consider the vector field: First, most of the neutron 
pairs inside a neutron star are in the triplet state. But second, since the graviton 
needs three degrees of freedom to become heavy, it would be more natural, if 
the symmetry breaking were done via one vector field instead of three scalar 
fields. 

The Lagrangian for a spin 1 field tp a which has both transversal and longi- 
tudinal degrees of freedom is given through Q : 

C = Q A^£> V + ^D^D^ - ~ ( M 2 - 2aR) ^ a - A (^ Q ) 2 ) » 

(47) 

where D a is the covariant derivative, 

A>V>/3 = - W (48) 

As before, we will redefine /j 2 — > fi 2 — 2aR to simplify the Lagrangian, 

£ = V9 {^D a ^ p D a ^ + °-D a ^D^ - ^ 2 M a - A (V Q <n 2 ) • (49) 

Again by variation with respect to the metric we obtain the following energy- 
momentum tensor 



T I 

kX 2 



g^D K ^D x ^ v + g^D^D^x - ~ 9kX g a0 g^ D a ^D^ 
+c (d k 4, x - i 9kX D fi i/A D a yj a - Q/i 2 + 2A^ a ) ^ x (50) 



i n , 



a 9k X I + XM a J M a - 

This looks a little complicated and a few simplifications are in order not to loose 
sight of what we are looking for. The first thing to notice is that the vacuum 
expectation value of the last term will always be zero, since above T c , \ip\ = 0, 
and below Tc, \ip\ 2 = —[i 2 /2\. Next we observe that the Christoffel symbols in 
the gravitational wave expansion p7| , 

= \l 1P (3aV + 9 PK P ~ d p K p ) + O (h 2 ) , (51) 

are proportional to d 1 h Xu , and thus we cannot expect them to give us a Meissner 
effect. Therefore, we will concentrate only on 



T I 

kA ~ 2 



g^d^Mu + g a0 d a ip K d^ x - i g KX g^g^d^d^ 



' [ d K ^ x - J g KX dp^A d a i> a - Q/i 2 + 2A^ Q V a ) > ; < '.v i ->2 1 
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We expand the metric = rf v — nh^ u and neglect terms in k 2 and re 3 . 
Furthermore assume that all the space-time dependance is in the phase, i.e., 

d K ^ = d K (\^\e i ^)=iip l ,d K ip. (53) 

Sorting with respect to rf v i\)^i\) v , we arrive at 

T K \ = 

j (((t?kA + Kh K \) r\ ati - K7] K xh afj ) d a (pdpip - 2d K ipdxip) jf v tp ll ip v 
+ Q (Kb?* - v af3 ) d a ipdpy - ~p 2 + 2A (Kh af> - v a0 ) $cSl>i\ Mx (54) 
( d K tpdxtp - \r] K xri a0 d a <pdp<p\ h^ u ^ u + c {KbT - rT) d„(pd K ^ v i(>x 



As in the scalar case, the above is going to be on the right-hand side of 
Einstein's equations ( f43| ) , and thus in the gravitational wave expansion 

d a d a h KX + (d K d x + Vnxd^) h a a = -8vrG T kA . 

Not unexpected, the energy-momentum tensor for the vector field is a little 
more complicated as the one for the scalar field. Direct interaction terms be- 
tween the graviton and the vector field like h^Tp^ip^ appear. However, the first 
line in equation (|54|) for the energy momentum tensor, i.e., 

t (({VkX + Kh K \) r] aP - n-q K \h al3 ) d a ipdpip - 2d K (pd\(p) if u ^^ip v , 

will show a phase transition. The vacuum expectation value of this term will 
be zero above the transition temperature and non-zero below, since 

IV'I 2 = for T > T c (55) 

= for T<T C . (56) 

Although a little more complicated than we may have wanted, the vector 
field too will give us a gravitational Meissner effect. This is very comforting in 
that it validates the assumptions in the previous sections. 

5.3 Discussion 

As a consequence of the Meissner effect, we expect the graviton to become 
massive. We can estimate the mass of the graviton inside a neutron star by 
taking the Compton wavelength formula 

m g =^— = 2.9 • l(T 47 kg = 1.6 ■ Kr n eV/c 2 , (57) 
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where Al is the London penetration depth and we used \l = 12km. As the 
graviton acquires a mass we then expect that the gravitational Coulomb poten- 
tial is to be replaced locally by a Yukawa type potential, with possible significant 
consequences for neutron star collapse, which are currently under investigation. 

6 Conclusion 

In this paper we tried to shed some light on the problem of breaking the gravi- 
tational symmetry. The approach taken is somewhat of a pedestrian in nature. 
As close a path as possible to observable objects and experiments was taken. 

The semiclassical London approach gives us the London penetration length. 
The number we obtain for neutron stars, \l = 12km, is intriguing and some- 
what unexpected, but is an indication that in neutron stars the breaking of the 
gravitational symmetry may play a significant role. 

The strength of the Ginzburg-Landau approach lies in that it is based on 
laboratory gravitational experiments. The COW experiment was the first one 
to show the direct influence of gravity on quantum mechanics. We go one step 
further and propose a minimal coupling scheme. From there everything else falls 
into place following in close analogy to the electromagnetic case. We predict 
again the Meissner effect, a gravitational Aharanov-Bohm effect, and for the 
triplet state a gravitational 'ferromagnetic' type phase. 

The last section addresses the problem from the gravitational wave expan- 
sion point of view. The spontaneous symmetry breaking character of the phe- 
nomenon becomes apparent. Becoming massive, the graviton will acquire three 
new degrees of freedom in the same way as the photon gets a new degree of 
freedom in the theory of superconductivity. 
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A Appendix 

A.l Remark on Units 

The treatment of the Ginzburg-Landau theory is significantly simplified if one 
follows the electromagnetic example and introduces the equivalent of Gaufiian 
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cgs-units. Similarly, as in the electromagnetic case, where the Coulomb force 

1 q 2 e 2 
F = * becomes F = 

47T£o r r 

where the two charges e and q are related through e = 4.80324 • 10~ 10 esu, one 
can introduce a new gravitational charge q, and 

2 1 
71% Cj 

F = G— y becomes F = -j, (58) 

where, for instance, 

<?at = x/Gtow = 1.35639 • 10~ 32 gsu (59) 

for a neutron with mAr = 1-66 • 10 _27 kg. 'gsu' stands for 'gravito-static units' 
and is gravitational equivalent of 'esu'. Table [l] lists a few units. 





units 


q = 


mass charge 


gsu = \/G kg 


q 2 . 


= (mass charge) 2 


(gsu) 2 = N m 2 


p 


density 


gsu / m 3 


j = 


current density 


gsu / m 2 s 


F = 


= force 


N = kg m / s 2 


g = 


electric field 


N / gsu 


H = 


= magnetic field 


N / gsu 


= 


potential 


N m / gsu 


A = 


= vector potential 


N m / gsu 



Table 1: Gravitostatic units. 



A. 2 Examples and Numbers 

Any gravitating body will be surrounded by a radial gravitoelectric field 

g = -vG T e r 



(60) 



where the bodies mass M is its gravitoelectric monopole moment. 

If the body is rotating, then in addition there will be a gravitomagnetic 
dipole field of the form 



H 



2VG 



(61) 



where the gravitomagnetic dipole moment is its spin angular momentum S 
Typical numbers for M and S are given in table §. 
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M 


S 


Earth 
Sun 

Neutron Star 


5.98- 10 2Y g 

1.99- 10 33 g 
2.8-10 33 g 


7-10 4U gem 2 / s 
1.7-10 48 gem 2 / s 
2-10 47 gem 2 / s 



Table 2: Some typical masses and angular momenti. 
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